Kondo impurities in nanoscopic systems: new confinement-induced regimes 
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We present results for Kondo impurities in nanoscopic systems. Using exact diagonalization in 
small clusters and Wilson's renormalization group we analyze an isolated system and a nanoscopic 
system weakly coupled to a macroscopic reservoir. In the latter case, new regimes not observed in 
macroscopic homogeneous systems are induced by the confinement of conduction electrons. These 
new confinement-induced regimes are very different depending on whether the Fermi energy coincides 
with the energy of a resonant state or lies between two quasi-bound states. 



PACS numbers: 72.I5.Qm, 73.22.-f 

The physics of nanoscopic and mesoscopic systems be- 
came of grate interest due to the possibility of making 
experiments in extremely small samples. Some exam- 
ples of nano and mesoscopic systems recently studied are 
smalL metallic and superconducting islands, El quantum 
dotsp nanotubesJ3 and quantum corralsD Depending on 
the type of system under study, different parameters can 
be controlled. 

The advances in nanotechnologies revived the interest 
in the Kondo effect, El one of the paradigms of strongly 
correlated systems. On one hand, Scanning Tunneling 
Microscopy (STM) allowed the direct measuremcri|t||ef lo- 
cal spectroscopic properties of Kondo impurities .raJ On 
the other hand it has been shown that quantum dots and 
single walled carbon nanotubes weakly coupled to con- 
tacts may behave as Kondo impurities generating new 
alternatives to study the phenomena.LJ 

In this work we report results for Kondo impurities in 
nanoscopic systems. Some experiments show that the ef- 
fects of confinement are important and recently Affleck 
and SimonEl proposed to use a mesoscopic system to mea- 
sure the Kondo screening length— Although the problem 
has been addressed theoreticallyJaO there are still many 
open questions.lia What is the temperature dependence 
of the Kondo screening in nanoscopic systems? What 
is the effect adding leads or coupling the system to a 
macroscopic reservoir? 

In the case of systems weakly coupled to a macroscopic 
reservoir, the one electron states acquire a finite lifetime 
and the spectrum consists of resonant states. This situ- 
ation is the most relevant to compare with experiments, 
however due to the intrinsic difficulties to formulate the 
problem is has not been analyzed in detail. As an exam- 
ple we may mention the case of quantum corrals, were 
the electrons are not well confined in the corral due to 
leaking of the wave function and hybridization with the 
bulk states. Consequently the number of electrons in 
the corral is not a good quantum number and the lo- 
cal density of states as a function of the energy presents 
well defined maxima-corresponding to the quasi-confined 
states of the corral£3 We show that in systems like the 
quantum corrals new regimes are induced by the confine- 
ment of conduction electrons. These new confinement- 



induced regimes are very different if the Fermi energy 
is at resonance or between two quasi-bound states. We 
also discuss the case of quantum dots in Aharonov-Bohm 
rings. 

We have performed numerical diagonalization in small 
clusters and implemented Wilson's numerical renormal- 
ization group (NRG)o Our starting point is the Ander- 
son model for magnetic impurities. The model Hamilto- 
nian reads :E3 

Ham = ^ £ d d a d <y + U d\d^d\di + ^ e v c\ a c va 

+ y]Vuci <T d <J + V*d\c va - fiiBS iz , (I) 



here the operator <ij. creates an electron with spin a at 
the impurity orbital with energy Ed and Coulomb re- 
pulsion U, cl a creates an electron in an extended state 



with quantum numbers v and a and energy e v . The 
last term represents the effect of an external magnetic 
field along the z-direction coupled to the impurity spin 
Si = S a T a . a id a i where r are the Pauli matrices. 

We start by analyzing the case of an isolated sys- 
tem. The properties of the system depend crucially on 
a few parameters, in particular we show below that in 
the Kondo limit a finite system presents quite a marked 
even-odd electron number asymmetry. 

We first consider the case of Hamiltonian (1) describ- 
ing a linear chain of N equivalent sites with the impu- 
rity at one end. In this case e v — —2tcos(vK/(N + I)), 
where t is the hopping matrix element in the chain and 
V v = Vo sin(wK/(N + I)) with v = 1, 2, . . . , N. We per- 
formed exact diagonalization in these finite systems using 
a Lanczos algorithm. 

For an even number of electrons in the system, the 
ground state is a singlet indicated as |0) and the expec- 
tation value of the impurity spin (0|iSyO) is zero. The 
impurity spin is screened by the free electrons and the 
zero temperature susceptibility is finite. We can esti- 
mate the characteristic energy scale fcsTfj- as the energy 
involved in the screening of the impurity spin. The mag- 
netization as a function of the external field B, shown in 
the inset of Fig.[l](a) , gives an estimation of Tk ■ We define 
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ksTx = HiB c where B c is the crossover field indicated in 
the inset. As t increases the mean density of conduction 
states decreases and consequently Tk decreases. As in- 
dicated in Fig.^(a) two regimes are obtained: IzbTk < A 
and UbTk > A where A ~ At/N is the level spacing. 
For the parameters used the crossover between these two 
regimes takes place at a hopping t cross ~ 0.27. 




FIG. 1. (a) Kondo temperature, as denned in the text, ver- 
sus hopping integral t for a chain of nine sites with eight par- 
ticles (U = 1.0, Vo = 0.3, and Ed = £f — 0.5). The dashed line 
is the level spacing A. Inset: impurity magnetization versus 
B. (b) Expectation value of the impurity spin versus t, for 
seven electrons, and the same parameters as in (a). 

For an odd number of electrons, the ground state is 
a Kramers spin-1/2 doublet, | ff) an d | JJ-)- The expec- 
tation value s = (-ft - \Si Z \ i[) = — (JJ- \Si Z \ i}-) is different 
from zero and, in the low temperature limit the impurity 
susceptibility diverges as \ = y^s 2 /ksT. In Fig.[j](b) the 
effective magnetic moment fj, = /XjS (we take /Ltj = 1) is 
plotted as function of t. For large A (large t) the im- 
purity spin is unscreened and s w 1/2 as expected for 
a gaped systemtS; for small A the impurity spin tends 
to be completely screened. The crossover between these 
two regimes takes place at t cross . 

Although these results were obtained for a very small 
system, they give the correct physical picture for zero 
temperature. While for A <C ksTx the impurity spin is 
essentially screened regardless of the parity of the elec- 
tron number, for A > ksTx an important even-odd 
asymmetry is obtained. In order to analyze larger sys- 
tems we have to resort to some approximate scheme. We 
have adapted the NRG to the present problem. The idea 
of the NRG is to obtain the quantum many-body states 
of the system on all energy scales, or length scales, in a 
sequence of steps. To do that, Wilson defined a basis of 
concentric wavefunctions in which the Hamiltonian takes 
the form of a linear chain. A truncated chain with n sites, 
described by an effective Hamiltonian H n , gives the cor- 
rect physics on an energy scale w n . The NRG transfor- 
mation T relates the Hamiltonians describing successive 
lower energy scales, H n+i = T[H n ], and leads to a sys- 
tematic way of calculating the thermodynamic properties 



at successive lower temperatures. 

Going back to our problem and for the sake of simplic- 
ity, consider a spherical metallic cluster of radius R c with 
the impurity at the center. The cluster is embedded in 
a bulk material with which it is weakly coupled through 
a large surface barrier. The Hamiltonian can be put in 
the form H = Ham + V(r, R c ) where the last term is 
a potential barrier placed at a distance R c from the im- 
purity. The Hamiltonian Ham is rewritten in Wilson's 
basis as schematically shown in Fig.||. The potential bar- 
rier is described by including a higher diagonal energy to 
the (n c + l)-shell centered around R c . Alternatively the 
barrier can be simulated with a smaller hopping matrix 
element in the region where V(r,R c ) is different from 
zero. We adopted this last description reducing the hop- 
pings by a factor a. 
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FIG. 2. Schematic representation of Wilson orbitals with 
the potential barrier and the corresponding linear chain. 

If the barrier is impenetrable (a — 0), the central clus- 
ter is isolated. The process of renormalization ends up 
after n c steps and with the obtained low energy states 
with a definite number of particles N e the partition func- 
tion Z is evaluated and the thermodynamic properties 
of the system are obtained. This procedure shows that 
the high temperature properties, evaluated with a num- 
ber of shells n < n c , the impurity behaves as in a bulk 
material. Only when the temperature is lowered and the 
discrete nature of the one-electron states becomes evi- 
dent the behavior of the system deviates from that of 
the macroscopic system. The impurity susceptibility can 
be evaluated for even and odd number of particles. As we 
discuss below, to describe even or odd number of parti- 
cles it is convenient to take even or odd-number of shells 
n c . The susceptibility is calculated as :liH 
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where the summation is done over the low energy states 
\v) with energies E v and P v — exp(-£^ /k B T) / 'Z '. The 
matrix elements {v\Si Z |£) have to be evaluated in a recur- 
sive way at each renormalization step. The susceptibil- 
ity reflects the thermodynamic properties of the system 
and we use it as an indicator of the degree of screening 
of the impurity spin. In Fig|j| the impurity susceptibil- 
ity is shown for different system sizes. For comparison 
the same quantity evaluated in the thermodynamic limit 
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(n c — > oo ) is shown. As n c decreases the characteristic 
energy separation between the one-electron states A in- 
creases. In the NRG this energy separation at the Fermi 
level is A ~ DA~ n "/ 2 where D is half the band width 
and A ~ 2 is Wilson's discretization parameter .E-3 
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FIG. 3. (a) Impurity magnetic moment squared for a finite 
system with an even number of particles as obtained with the 
NRG, the parameters are D = 1, U = 0.02, Ed = —0.01, and 
V = 0.01. The thick line corresponds to the thermodynamic 
limit, thin lines, from right to left, to n c = 46, 50, 56, and 
60. (b) Same as in a for (a) system with an odd number 
of particles, the thin lines from top to bottom correspond 
to n c = 45, 51, 55, 61, 65, and 71. Inset: Zero temperature 
impurity magnetic moment squared versus log(i? c /£|?). 

In Figl|(a) the effective impurity magnetic moment 
squaredEfT/U 2 = k B T\ for an even number of electrons is 
shown as a function of the temperature. When the tem- 
perature is of the order of A the behavior deviates from 
that of the macroscopic sample and /i 2 decreases expo- 
nentially due to the gaps induced by confinement. For 
an odd number of electrons the behavior is illustrated in 
Fig.|(b). In this UbT approaches A the screen- 

ing of the impurity spin ends up and a finite magnetic 
moment survives down to zero temperature, the ground 
state is a spin- 1/2 doublet. In the inset the low tem- 
perature magnetic moment versus log(i? c /£|?) is shown 
with £|? being the Kondo screening length (see below), 
the behavior shows the same tendency of that of Fig.|I](b) 
corresponding to exact results of the impurity in a linear 
chain. Here the crossover between an almost unscreened 
impurity spin to an almost completely screened spin oc- 
curs for A ~ k B T™. 

The behavior of the system with a finite barrier is 
richer and more interesting. Now the central cluster is in 
contact with a reservoir and the number of electrons in 
the cluster is not a good quantum number. In the numer- 
ical NRG approach the Fermi energy is set to zero. The 
one particle spectrum in the absence of the impurity has 
a state at zero energy for even iterations n and two states 
at approximately ±DA - "/ 2 for odd n. Then by taking 
the barrier enclosing an even or odd number of shells, n c , 
we describe a system with the Fermi level at the energy 



of a resonant state (at resonance) or between two of them 
(off resonance) respectively. These two situations corre- 
spond approximately to an even and odd average number 
of electrons in the cluster respectively. 

In Fig.|(a) the impurity magnetic moment /i 2 is shown 
for the a system at resonance. Only when the temper- 
ature is lowered and reaches the value fc^T « A the 
structure in the local density of states becomes impor- 
tant. At low temperatures the universal Kondo behavior 
is recovered with a Kondo temperature T|? higher than 
Tp?, corresponding to a higher density of states. The 
crossover between the high and low temperature regimes 
occurs in a small temperature interval leading to a rapid 
decrease of the magnetic moment \j? . 

For the system off resonance, as fcgT approaches A the 
magnetic moment saturates leading to a plateau in the 
curve /i 2 versus T as shown in Fig.^(b). Only at lower 
temperatures the universal Kondo behavior is recovered 
with a Kondo temperature T^ that is lower than Tg?. 
The plateau corresponds to a new regime of partially 
screened magnetic moment. In this regime, as the tem- 
perature is lowered, the thermodynamic properties do 
not change and are similar to those of a finite system 
with an even number of particles. The central cluster 
has a magnetic moment and it is the whole cluster that 
behaves as a Kondo impurity or a "Kondo grain" when 
the temperature is lowered. For a given potential bar- 
rier height the size of the grain can be changed: systems 
with small R c , a large A, present a wide plateau with 
a partially-screened magnetic moment /i 2 s ~ 1/4. As 
R c increases, both the width ST of the plateau where the 
partially-screened magnetic moment fi ps is stable and the 
value of iip S decrease. For a given size, the ST increases 
with the barrier height as illustrated in Fig.||(c). 

In summary, for fc^T^ ^> A the fine structure (on 
the scale of TgP) of the density of states does not change 
the properties of the system and the universal Kondo be- 
havior is recovered. For ksT^ <C A the local density 
of states, on the scale of T2?, is a smooth function and 
the universal Kondo behavior is obtained at low temper- 
atures. In this case the system at resonance has a Kondo 
temperature T R much larger than the corresponding to 
the system out of resonance T^ R , however in both cases 
the high temperature susceptibility is non universal. For 
hsT^? ~ A new confinement induced regimes are ob- 
served: for the system at resonance, as the temperature is 
lowered, there is a rapid decrease in the magnetic moment 
/j, 2 when fcfiT ~ A; for the system off resonance there is 
a new partially-screened regime that leads to a plateau in 
the temperature dependence of fi 2 . This plateau can be 
interpreted as the free-spin regime of the whole Kondo 
grain that has magnetic moment partially localized at the 
impurity. At very low temperatures the electrons that 
are outside the cluster complete the screening. The tem- 
perature at which this occurs depends on the coupling 
between nanoscopic system and the bath. The condition 
for the existence of these new regimes k B T^ ~ A can be 
put in the form £|? ~ R c where = Hvf /ksT^ is the 
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Kondo screening length. 
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new effects are to be taken into account. In this case 
the quantum dot behaves quite differently if the Fermi 
level lies at resonance or off resonance, corresponding to 
a maximum or a minimum in the transmittance respec- 
tively. The magnetic field can continuously shift the res- 
onances and if the system is in the temperature interval 
, even at constant temperature the QD 
continuously evolves from a low temperature regime (T 
<T^) to a high temperature regime (T^ R <T ). More- 
over, the energy of the resonances depend not only on 
the external field but also on the phase shift introduced 
by the QD that is sclfconsistently adjusted. On the light 
of our analysis, in this regime the AB oscillations can not 
be interpreted directly in terms of a Kondo impurity at 
a fixed temperature ratio T/Tk- 
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FIG. 4. Impurity magnetic moment squared for a sys- 
tem with a finite potential barrier and the Fermi energy at 
resonance (a) and off resonance (b) and (c). Parameters are 
the same as in Fig. (3). In (a) and (b), a — 0.05 and dif- 
ferent cluster sizes are shown (n c ~ 50, 62, 72 in (a) and 
n c — 49, 61, 71 in (b)). In (c) the potential barrier is changed 
for a fixed size: n c = 61, a = 0, 0.05, 0.1, 0.3, and 1, (from 
left to right). In (b) and (c) the thin lines that extrapolate 
to a finite value correspond to the isolated system. In a) and 
b) the insets show a typical local density of states of the con- 
duction electrons at the impurity site. 

For the case of Co impurities in quantum corrals built 
on the Cu (111) surface, taking T£? = 50K and using 
the parameters of the Cu (111) surface band we estimate 
£|? ~ 400A which is larger than the typical corral ra- 
dius. Consequently if the Fermi level is out of resonance 
we expect the Kondo screening to occur only at very low 
temperatures. Larger corrals could be built or the ef- 
fective barrier height decreased by removing some atoms 
from the corral fence and in this way control the effective 
Kondo temperature of the impurity. 

Our results may have important consequences for the 
interpretation of .experiments with QD in Aharonov- 
Bohm (AB) ringsBQ In the AB configuration, the local 
density of states inside the ring has the structure of res- 
onant states. A direct consequence of this structure are 
the oscillations in the transmittance of a ring as a mag- 
netic field threading the circuit is varied. For rings built 
on GaAs as in some of the experimental setups the energy 
separation between resonances A could be larger than 
one Kelvin, in fact at IK oscillations are still observed 
indicating that at this temperature thermodynamic fluc- 
tuations do not wash out the density of states structure. 
If the Kondo temperature of a quantum dot inserted in 
one of the ring arms is of the order or smaller than A, 
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